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Full-sky maps of the cosmic microwave background temperature reveal a 7% asymmetry of fluctu¬ 
ation power between two halves of the sky. A common phenomenological model for this asymmetry 
is an overall dipole modulation of statistically isotropic fluctuations, which produces particular off- 
diagonal correlations between multipole coefficients. We compute these correlations and construct 
corresponding estimators for the amplitude and direction of the dipole modulation. Applying these 
estimators to various cut-sky temperature maps from Planck and WMAP data shows consistency 
with a dipole modulation, differing from a null signal at 2.5cr, with an amplitude and direction con¬ 
sistent with previous fits based on the temperature fluctuation power. The signal is scale dependent, 
with a statistically significant amplitude at angular scales larger than 2 degrees. Future measure¬ 
ments of microwave background polarization and gravitational lensing can increase the significance 
of the signal. If the signal is not a statistical fluke in an isotropic Universe, it requires new physics 
beyond the standard model of cosmology. 

PACS numbers: 98.70.Vc, 98.80.Es, 98.80.Jk 


I. INTRODUCTION 

The statistical isotropy of the cosmic microwave background (CMB) at large angular scales has been questioned 
since the first data release of the WMAP satellite [T]. Independent studies performed on different WMAP data 
releases mn show that the microwave temperature sky possesses a hemispherical power asymmetry, exhibiting more 
large-scale power in one half of the sky than the other. Recently, this finding has been confirmed with a significance 
greater than Scr with CMB temperature maps from the first data release of the Planck experiment The power 
asymmetry has been detected using multiple techniques, including spatial variation of the temperature power spectrum 
for multipoles up to 1 = 600 [H] and measurements of the local variance of the CMB temperature map [3 E] ■ For 
I > 600, the amplitude of the power asymmetry drops quickly with I [HE]. 

A phenomenological model for the hemispherical power asymmetry is a statistically isotropic sky 0(n) times a 
dipole modulation of the temperature anisotropy amplitude: 

0(n) = (1 + n-A)0(n), (1) 

where the vector A gives the dipole amplitude and sky direction of the asymmetry m- This phenomenological model 
has been tested on large scales (Z < 100) with both WMAP [TTJ [T^] and Planck (0, hereafter PLK13) data, showing 
a dipole modulation with the amplitude |A| ~ 0.07 along the direction {£,b) ~ (220°,—20°) in galactic coordinates, 
with a significance at a level > 3cr. Further analysis at intermediate scales 100 < I < 600 shows that the amplitude 
of the dipole modulation is also scale dependent [T5] . 

If a dipole modulation in the form of Eq. Q is present, it induces off-diagonal correlations between multipole 
components with differing I values. Similar techniques have been employed to study both the dipole modulation 
[SHISl and the local peculiar velocity [nMO]. In this work, we exploit these correlations to construct estimators 
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for the Cartesian components of the vector A as function of the multipole. These estimators are then applied to 
publicly available foreground-cleaned Planck CMB temperature maps. We constrain the scale dependence over a 
multipole range 2 < I < 600, as well as determine the statistical significance of the observed geometrical conhguration 
as function of the multipole. Throughout this analysis, we adopt realistic masking of the galactic contamination. We 
test our findings against possible instrumental systematics and residual foregrounds. 

In the following Section, we derive estimators for the dipole modulation components and their variances for a 
cosmic-variance limited CMB temperature map. Section [ITT] presents and tests a pipeline for deriving these estimators 
from observed maps, showing how to correct for partial sky coverage. Using simulated CMB maps, we estimate the 


covariance matrix of the components of the dipole vector, as well as testing for possible systematic effects. Section IV 


describes the Planck temperature data we use to obtain the results in Sec. |V] We estimate the components of the 
dipole modulation vector and assess their statistical signihcance, finding departures from zero at the 2 to 3 u level. 
The best-fit dipole modulation signal is an unexpectedly good ht to the data, suggesting that we have neglected 
additional correlations in modeling the temperature sky. We also perform a Monte Carlo analysis to estimate how 
the dipole modulation depends on angular scale, confirming previous work showing the power modulation becoming 
undetectable for angular scales less than 0.4°. Finally, Sec. |VI| gives a discussion of the significance of the results and 
possible implications for models of primordial perturbations. 


II. DIPOLE-MODULATION-INDUCED CORRELATIONS AND ESTIMATORS 


Assuming the phenomenological model described by Eq. 0. the dipole dependence on direction can be expressed 
in terms of the I = I spherical harmonics as 

n • A = 2Y^(A+Yi_i(ri) — A_Yi+i(n) -|- A^yio(n)) (2) 

with the abbreviation A± = {A^^iAy) / ^/2. Expanding the temperature distributions in the usual spherical harmonics, 

0(n) = ^ aimV/m(n), 0(n) = ^ aimYim{^), (3) 

Im Im 

with the usual isotropic expectation values 


{o*lmO‘l'm') = Cl5u'5ram'■ (4) 

The coefficients must satisfy = (—and a*^ = {—l)"^ai-m because the temperature held is real. The 
asymmetric multipoles can be expressed in terms of the symmetric multipoles as 


aim = azm + ^ J dhYi_min)Yi:m'{n) [A+Yi_i(ri) - A_Yi+i(n) -H A2Yio(n)]. 

I'm' 

The integrals can be performed in terms of the Wigner 3j symbols using the usual Gaunt formula, 

f _ . /(2^1 + 1)(2/2 + 1)(2Z3 + 1) ( h h h \ fh h h 

y (n) (n) - y waj (^0 0 0 


( 5 ) 


( 6 ) 


Because h = 1 for each term in Eq. ([^ , the triangle inequalities obeyed by the 3j symbols give that the only nonzero 
terms in Eq. ([^ are I' = I or I' = I ± 1. For these simple cases, the 3j symbols can be evaluated explicitly. Then it is 
straightforward to derive 


{ai+lm±iaim) — ^Ui) 


{l±m + 2){l ± m -I- 1) 
(2/ + 3)(2Z + l) ’ 


mdlm) — {Cl + Cl+i) 


{I — m + 1){1 + m + 1) 
{21 -|- 3)(2^ -|- 1) 


( 7 ) 

( 8 ) 


These off-diagonal correlations between multipole coefficients with different I values are zero for an isotropic sky. This 
result was previously found by El: and represents a special case of the bipolar spherical harmonic formalism |2I| . 
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It is now simple to construct estimators for the components of A from products of multipole coefficients in a map. 


\/2ReA+ and Ay = •\/2ImA_|_, we obtain the following estimators: 

111’’" ^ 

1 {21 3)(2/ + 1) 

(R,6 777,_|_iR,G QjI 

"J'- - Cl + Ci+, \ 

1 {l + m + 2){l + m+l) 


1 [21 3)(2/ + 1) 

(Re fo+i,„+ilm ai 

Vhm - Cl + Ci+I \ 

I {1 m 2){l + m + 1) 

14 1''' ^ 

1 [21 + 3)(2/ + 1) 

(RG 777,R,G 


/ {1 + m + 1){1 — m -1- 1) 


(9) 

( 10 ) 

( 11 ) 


where ajm’s are calculated from a given (real or simulated) map and C;’s are estimated from the harmonic coefficients 
of the isotropic map Ci = {21 + 1)~ ^ We argue that for small values of the dipole vector A and (more 

important) for a nearly full-sky map ^ t X) lo/mP- This assumption has been tested for the kinematic dipole 

modulation induced in the CMB due to our proper motion, showing that the bias on the estimated power spectrum 
is much smaller than the cosmic variance error for nearly full-sky surveys [20] , Such estimators, derived under the 
constraint of constant dipole modulation, can be safely used for the general case of a scale-dependent dipole vector 
A by assuming that A.{1) ~ A(/ -|- 1). This requirement is trivially satisfied by a small and monotonically decreasing 
function A{1). 

To compute the variance of these estimators, assume a full-sky microwave background map which is dominated by 
cosmic variance; the Planck maps are a good approximation to this ideal. Then aim is a Gaussian random variable 
with variance af = Ci. The real and imaginary parts are also each Gaussian distributed, with a variance half as large. 
The product x = Redi+im-i-iR-6®/m, for example, will then have a product normal distribution with probability 
density 


P{x) = -Ao ^ 

TrCTifTi+i \(T;cr/+i 




( 12 ) 


with variance cr^ = where Kq{x) is a modified Bessel function. By the central limit theorem, a sum of 

random variables with different variances will tend to a normal distribution with variance equal to the sum of the 
variances of the random variables; in practice the sum of two random variables each with a product normal distribution 


will be very close to normally distributed, as can be verified numerically from Eq. (12). Therefore we can treat the 
sums of pairs of aim values in Eqs. (|^ „2 


errors for the estimators as 


to (|ll|) as normal variables with variance (j'faf^-^/2, and obtain the standard 


. 1 — 1 ~ / (2/-I-3)(2Z-I-1) 

^\im y vhm-2{lJ^rn + 2){l + m + iy 

^1 / (2/ + 3)(2Z + l) 

1 2\j 2{l + m + l){l-m + iy 


(13) 

(14) 


with the approximation ~ Ci. 

Eor a sky map with cosmic variance, each estimator of the components of A for a given value of I and m will have a 
low signal-to-noise ratio. Averaging the estimators with inverse variance weighting will give the highest signal-to-noise 
ratio. Consider such an estimator for a component of A which average all of the multipole moments between I = a 
and I = b\ 


l—a m— — l 

l—a m— — l 

[•4.1='’•EE 
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which have standard errors of 



' b 1 

-1/2 


Cj. (Ty - 

.l—a m——l 

= 

■ 2 

— (6 + 0 + 2) (6 — o + 1) 
o 



.l—a m—0 

-1/2 

4(6 — 0 + 1) [o(26 + 3)(o + 6 + 4) + (6 + 2)(6 + 3)] 

az — 


3(2a-kl)(26-p3) 



(18) 

(19) 


The sum over m for the z estimator and error runs from 0 instead of —I because [Az\i-m = [Az]irm but the values are 
distinct for the x and y estimators. 

While the Cartesian components are real Gaussian random variables, such that for isotropic models ([^x]) = 
{[Ay\) = {[A^]) = 0, the amplitude of A is not Gaussian distributed. Instead, it is described by a chi-square 
distribution with 3 degrees of freedom, which implies (|Ap) ^ 0 and p(|Ap = 0) = 0, even for an isotropic sky. For 
this reason, we consider the properties of the dipole vector A as function of the multipole, considering each Cartesian 
component separately. 


III. SIMULATIONS AND ANALYSIS PIPELINE 


The estimators in Eq. ([^-(11 1 are clearly unbiased for the case of full-sky CMB map. However, residual foreground 
contaminations along the galactic plane as well as point sources may cause a spurious dipole modulation signal, which 
can be interpreted as cosmological. Such highly contaminated regions can be masked out, at the cost of breaking the 
statistical isotropy of the CMB held and inducing off-diagonal correlations between different modes. The effect of the 
mask, which has a known structure, can be characterized and removed. 


A. Characterization of the Mask 


For a masked sky the original aim are replaced with their masked counterparts: 


aim = J dne{n)W{n)Yi* 


( 20 ) 


where W{h) is the mask, with 0 < kF(n) < 1. In this case, Eq. Q does not hold, meaning that even for a statistical 
isotropic but masked sky the estimators in Eq. (|^-([TT|) will have an expectation value different from zero. This 
constitutes a bias factor in our estimation of the dipole modulation. 

If we expand Eqs. using the dehnition of the harmonic coefficients in Eq. ([^ , it is clear that if a primordial 

dipole modulation is present, the mask transfers power between different Cartesian components. Under the previous 
assumption A.{1) ~ A(/ -|-1), the Cartesian components i,j = x,y, z of the dipole vector can be written as 


— Aji imAi i -f Mj im 


( 21 ) 


where [Aj]im is the estimated dipole vector for the masked map, and Aji^im and Mjjm are Gaussian random numbers 
determined by the d/m, so dependent only on the geometry of the mask. For unmasked skies, these two quantities 
satisfy (Aj//m) = and (Mj/m) = 0, ensuring that the expectation value of our estimator converges to the true 
value. 


Using Eq. (211, we can define a transformation to recover the true binned dipole vector from a masked map: 


[A,] = A-/([A,] - M,) 


( 22 ) 


where [Aj\ is the binned dipole vector estimated from a map, and A^/ and Mj are the expectation values of Aji^im 
and Alj im, binned using the prescription in Eqs. (15|-(17). For each Cartesian component we divide the multipole 
range in 19 bins with uneven spacing, Al — 10 for 2 < [< 100, Al — 100 for 101 < I < 1000. For a given mask, the 
matrix Aji and the vector Mj can be computed by using simulations of isotropic masked skies. We use an ensemble 
of 2000 simulations, and we adopt the apodized Plauick U73 mask, following the procedure adopted by PLK13 for the 
hemispherical power asymmetry analysis. For the rest of this work, all estimates of the dipole vector are corrected 
for the effect of the mask using Eq. (22). 
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Cartesian components 


Cartesian components 


FIG. 1. Correlation matrices for the Cartesian components of the dipole vector. These matrices are estimated using 2000 random 
simulated skies masked with the apodized Planck U73 mask. The ordering of the components follows the convention defined for 
the dipole vector. (Left panel) Isotropic skies (A=0). (Right panel) Difference between the correlation matrices for modulated 
skies, generated using a constant dipole vector across multipoles of magnitude |A| = 0.07 and direction {i,b) = (220°,—20°), 
and isotropic skies. 


B. Simulated Skies 

We generate 2000 random masked skies for both isotropic and dipole modulated cases. For the latter, we assume an 
scale-independent model with amplitude |A| = 0.07, along the direction in galactic coordinates (/,6) = (220°,—20°). 
We adopt a resolution corresponding to the HEALPi3|^ [22] parameter NSIDE = 2048, and we include a Gaussian 
smoothing of FWHM = 5' to match the resolution of the available maps. The harmonic coefficients dim are then rescaled 
by \/5/, where the power spectrum is calculated directly from the masked map. These normalized coefficients (for 
both isotropic and dipole modulated cases) are then used to estimate the components of the dipole vector. 

These simulations also serve the purpose of estimating the covariance matrix C. From Eqs. we expect 

different Cartesian components to be nearly uncorrelated, even for models with a non-zero dipole modulation, for 
full-sky maps. We confirm this numerically with simulations of unmasked skies. For masked skies. Fig, ([^ shows the 
covariance matrices. The left panel shows the case for isotropic skies with no dipole modulation. The presence of the 
mask induces correlations between multipole bins at scales 100 < I < 500, and also between the largest scales I < 40 
with all the other multipole bins. However, because of the apodization applied to the mask, the correlation between 
bins never exceeds 25%. For comparison, we also show the difference between the correlation matrices for the case 
of dipole modulated and isotropic skies (right panel). This is consistent with random noise, which demonstrates that 
the covariance matrix does not depend significantly on the amplitude of the dipole modulation. 


C. Bias Estimates 


We determine the mean bias in reconstructing the dipole modulation vector A from a masked sky by computing 
the mean value of all three Cartesian components reconstructed from 2000 simulations, for both isotropic and dipole- 
modulated skies. In both cases, the residual bias vector has components Ai > 0, with an amplitude of the first bin 
of each Cartesian component below 6 x 10“^. For the isotropic case, the amplitude of the bias is strongly decreasing 
with multipole {\A{1 = 60)| = 3.8 x 10“^), corresponding to 0.5% to 2% of the cosmic variance error for the entire 
multipole range considered. Therefore, the analysis procedure on masked skies does not introduce a statistically any 
significant signal which could be mistaken for dipole modulation. 


In the case of dipole-modulated simulations with dipole amplitude A = 0.07 consistent with PLK13 (see Sec. Ill C I, 


the amplitude of the bias for each Cartesian component is a constant for all multipoles. This indicates that the bias 
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FIG. 2. Measured Cartesian components of the dipole vector from the SMICA Planck map as a function of the central bin 
multipole /center- The amplitudes are multiplied by y/l to enhance visibility of the signal at higher multipoles. The Icr errors 
are the square roots of the covariance matrix diagonal elements. Data at Z > 600 is not used in our statistical analyses. 


follows the underlying model, and the determination of the scale dependence of the true dipole vector will not be 
affected by such a bias. For this specific case, the amplitude is always < O.Scr when compared to the cosmic variance, 
specifically < O.lcr for I < 100. However, this simulated case is unrealistic. We do not expect such a big amplitude 
for the dipole vector at small scales, so the simulated case overestimates the actual bias. 


IV. MICROWAVE SKY DATA 


We consider a suite of six different foreground-cleaned microwave background temperature maps|^ SMICA, NILC 
COMMANDER-Ruler H and SEVEM from the first Planck data release [23], and two others processed with the LGMCI0 
component separation technique by |2l]. The LGMCA-PRl and LGMCA-WPRl are based on Planck data only and 
Planck-I-WMAP9 data respectively, allowing a non trivial consistency test between these two experiments. Each of 
these maps uses a somewhat different method for separating the microwave background component from foreground 
emission, allowing us to quantify any dependence on the component separation procedure. 

Asymmetric beams and inhomogeneous noise may create a systematic dipolar modulation in the sky. In order 
to test this possibility, we analyze the 100 publicly available FFP6 single-frequency simulated maps released by the 
Planck team. Specifically, we process the simulations for the 100, 143 and 217 GHz channels with our analysis code. 
The maximum likelihood analysis shows a bias on small scales, although the values are always less than 0.6 times the 
cosmic variance for each multipole bin. Considering only the first 15 bins (Zmax = 600) gives a result consistent with 
the isotropic case, with a p-value larger than 0.1. The source of the small-scale bias is not yet known, but we simply 
ignore multipoles I > 600 in the present analysis of Planck data. 


^ http://wiki.cosmos.esa.int/planckpla/index.php/Main_Page 
® http: //WWW. cosmostat.org/product/lgmca_cmb/ 
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Multi pole 

FIG. 3. Measured amplitude of the dipole vector from the SMICA Planck map. The black dashed line shows the model for the 
statistical isotropic case A=0. The Icr errors are estimated from the 16*^ and 84*^ percentiles of the distribution of the dipole 
vector amplitudes, calculated from sky simulations processed the same way as the data. 


V. RESULTS 

Figure ([^ shows the measured values of the Cartesian components of the dipole vector, using the SMICA map. 
Similar results are found for the other foreground-cleaned maps, and a direct comparison is shown in Sec. |V A| 
Figure ([^ displays the amplitude of the dipole vector compared with the mean value (black dashed line) obtained 
from isotropic simulations; as pointed out in Sec. [l^ the expectation value of the amplitude of the dipole vector is 
different from zero even for the isotropic case. 

The data clearly shows two important features: 

• The amplitudes of the components of the estimated dipole vector are decreasing functions of the multipole 1. 

• The X and y components have a negative sign, which persists over a wide range of multipoles; the z component 
is consistent with zero. This indicates that the vector is pointing in a sky region (180° < i < 270°, b ~ 0), in 
agreement with previous analyses. 

We further characterize these basic results in the remainder of this Section. 


A. Geometrical Test 


First, we test how likely the observed geometrical configuration of the dipole vector is in an isotropic universe. To 
achieve this goal, we need to define a quantity which preserves the information on the direction of the dipole vector 
(i.e. statistics linear in the variables Ai). In addition, the Cartesian components have to be weighted by the cosmic 
variance, ensuring that our statistics is not dominated by the first bins. Therefore, we define the following quantity 


3N 


a 




2=1 


IIIB 


(23) 


and [4l]j=i,...,3Af 


where (C“^)ij are the components of the inverse of the covariance matrix calculated in Sec. 
are the three Cartesian components of the binned dipole vector (up to the TV**' bin) estimated either from a simulated 
map of measured data. For an isotropic universe, we expect the three Cartesian components to sum up to zero, 
such that (a) = 0 for any choice of /max- This will not be the case if the underlying model is not isotropic (i.e. 
the expectation values of the Cartesian components are different from zero). In Fig. we plot the values of the a 
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FIG. 4. The a—parameter from Eq. (231, scaled by the standard deviation cr{a), as function of the maximum multipole 
considered /max- The colored solid lines are the results from CMB data, showing remarkable consistency between different 
foreground cleaning methods. (Left panel) The shaded bands are estimated using simulations of isotropic masked skies. The 
distribution a—parameter is Gaussian with (a) = 0. (Right panel) The shaded bands are estimated using simulations of dipole 
modulated masked skies. The dipole modulation model is ^ = 0.07, {i,b) = (220°,—20°). The confidence regions (colored 
band) are estimated using percentiles, such that Tier = [15.87*'*, 84.13*'"], ±2 ct = [2.28"'*, 97.72*'*] and ±3cr = [0.13“*, 99.87*'*]. 


parameter as function of the maximum multipole considered in the analysis Imaxi rescaled by the standard deviation 
(j{a) determined from the simulations of isotropic skies. The left panel shows the comparison between the CMB data 
for all 6 foreground-cleaned maps, and the simulations for the isotropic case. The measured rescaled a parameter has 
a value that is discrepant from a = 0 at a level of 2cr < a < Scr. This discrepancy is maximized for I 70 — 80, which 
corresponds to what has been previously probed by PLK13. 

The right panel of Fig. compares the measured signal with simulations of dipole modulated skies, using the 
covariance matrix C calculated from the anisotropic simulations. This test confirms that the signal averaged over 
multipoles ^ 70 is consistent with the model proposed by PLK13 (assumed in our anisotropic simulations). However, 
the results are not consistent with a scale-independent dipole modulation, and the amplitude of the dipole modulation 
vector must be strongly suppressed at higher multipoles. 


B. Model Fitting 


Consider a simple power-law model for the dipole modulation defined by 4 parameters: 


=4( —I cos&cos£, 

V60/ 

(24) 

/ 1 \” 


=4 — cos&sin£, 

V60/ 

(25) 

/ 1 


= ->l(^) smi, 

(26) 


where A is the amplitude of the dipole vector at the pivot scale of Z = 60, n is the spectral index of the power law, b 
is the galactic latitude, and i is the galactic longitude. We use a Gaussian likelihood £, such that 


In/: = = -i([4l.] - [A,rfiC-%{[A,] - [A,r) 


(27) 


where [Aj\ are the estimated components from the Planck SMICA map, are the components of the assumed 

model properly binned using Eqs. ©-([Tg, and C is the covariance matrix for a dipole-modulated sky displayed in 
Fig. 0 The parameter space is explored using the Markov chain Monte Carlo sampler emce^ [^, assuming flat 
priors over the ranges {A,n} = {[0,1], [—2,2]}. Tablejl], displays the results for different thresholds of Zmax- 
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TABLE I. Best-fit values of the amplitude A, spectral index n and direction angles {£, b) for the dipole vector, as function of the 
maximum multipole Zmax- The best fit values are corresponds to the 50^^ percentile of the posterior distribution marginalized 
over the other parameters. The errors corresponds to the 16*^ and 84*^ percentiles. For the first case, we consider a model 
with spectral index n = 0 and /max = 60, which can be compared with PLK13 findings. Values of the corresponding to the 
best fit model, as well as to the isotropic case, are also displayed with corresponding number of degrees of freedom i/. 


^max 

A 

n 

fori 

41 

xLinH 

xi=o(+ 

60 

U.UDO o.o30 

— 

-IOI 2 I 

+218tlt 

10.3 (15) 

19.5 (18) 

200 

U.UO4_0.0l6 

-0.54«:3| 

y+16 
“ ‘ -16 

+21+1® 

16.4 (29) 

30.8 (33) 

300 

^•^^^-0.014 

-o.68«:fg 

_11 + 18 
-1-1-16 

+21ltfo 

18.1 (32) 

31.4 (36) 

400 

U.UZ/ _Q Q;,^4 

-0.74«;fg 

_Q+19 

1^-18 

+212^11 

19.3 (35) 

31.9 (39) 

500 

0 0S1 

U.UOi_o.oi3 

-O.Oltg.fs 

-4+13 

^-13 

+20711^ 

22.5 (38) 

40.2 (42) 

600 

0 031+0 011 
U.U01 Q Q;L2 

-0.64t°:« 

1+13 

-1-14 

+20911® 

23.9 (41) 

41.9 (45) 


In the restricted case considering only low multipoles I < 60 and a flat spectrum n = 0, our best-fit model agrees 
at the 1(7 level with previous analysis by PLK13, for both amplitude and direction. 

If n is allowed to vary, the amplitude A of the dipole vector at the pivot scale of 60, as well as the spectral index n, 
are perfectly consistent for three different Zmax thresholds. The amplitude is different from the isotropic case A = 0 at 
a level of 2a, and the scale-invariant case n = 0 with = 400 is excluded at greater than 3a significance. The value 
of the galactic longitude £ is stable to a very high degree, whereas the value of galactic latitude b indicates (although 
not statistically significant) a migration of the pointing from the southern hemisphere to the northern one. This is 
expected due to the effect of the kinematic dipole modulation induced by the proper motion of the solar system with 
respect to the microwave background rest frame [niEe]. This effect has been detected by Planck m, and results in 
a dipole modulation in the direction {£,b) = (264°, 48°) detectable at high 1. 

The dipole model is a better fit to the data than isotropic models. Both the Aikake Information Criterion (AIC) 
and the Bayes Information Criteron (BIG) [SJ show sufficient improvement in the fit to justify the addition of four 
extra parameters in the model. In the specific case of the AIC the dipole model is always favored. The improvement is 
calculated by the relative likelihood of the isotropic model with respect to the dipole modulated case. This is defined 
as exp((AICmin — AICa=o)/2), where the AIC factor is corrected for the finite sample size, and it corresponds to 
0.48,0.083,0.13,0.18,0.013 and 0.011 for the models considered in Table|lj In the case of the BIC, the corresponding 
values of BICmin — BICa=o = —0.5, —0.4,1.0, 2.1, —2.8 and —2.8. The BIC indicates that the dipole modulation is 
favored only for the cases with £max > 400, where the parameters are better constrained. According to [55], the 
improvement even though positive is not strong because —6 < BICmin — BICa=o < —2. 

For the dipole-modulated model, the value of is substantially lower than the degrees of freedom. This suggests 
that either the error bars are overestimated, or the data points have correlations which have not been accounted for 
in the simple dipole model. Since the errors are mostly due to cosmic variance on the scales of interest, the error 
bars cannot have been significantly overestimated. Therefore, our results may point to additional correlations in the 
microwave temperature pattern beyond those induced by a simple dipole modulation of Gaussian random anisotropies. 
The correlations are unlikely to be due to foregrounds, since the results show little dependence on different foreground 
removal techniques. 


VI. CONCLUSIONS 

The microwave sky seems to exhibit a departure from statistical isotropy, due to half the sky having slightly more 
temperature fluctuation power than the other half. This work shows that the temperature anisotropies are consistent 
with a dipolar amplitude modulation, which induces correlations with multipole coefficients with I values differing by 
one. At angular scales of a few degrees and above, the correlations define a dipole direction which corresponds to the 
orientation of the previously known hemispherical power asymmetry, while at smaller scales the direction migrates 
to that of the kinematic dipole. Our results show that a dipole modulation is phenomenologically a good description 
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of the power asymmetry, but that the modulation must be scale dependent, becoming negligible compared to the 
kinematic dipole correlations [I71II9] on angular scales well below a degree. 

The statistical significance of these multipole correlations is between 2(j and 3cr compared to an isotropic sky, 
with the error dominated by cosmic variance. The maximum signal appears at scales Z ;< 70 as seen previously by 
PLK13. We also find an unusually low scatter in the dipole component estimates as a function of scale, given the 
cosmic variance of an unmodulated Gaussian random field, suggesting that the microwave temperature sky may have 
additional correlations not captured by this simple model. 

On the largest scales of the universe, simple models of Inflation predict that the amplitude of any dipole modulation 
due to random perturbations in a statistically isotropic universe should be substantially smaller than that observed. 
This departure from statistical isotropy may require new physics in the early Universe. One possible mechanism 
is a long-wavelength mode of an additional field which couples to the field generating perturbations P5H57] . If the 
mode has a wavelength longer than the current Hubble length, an observer sees its effect as a gradient. The field 
gradient modulates background physical quantities such as the effective inflaton potential or its slow-roll velocity. The 
required coupling between long and short wavelength modes can be accomplished in the context of squeezed-state 
non-Gaussianity |49fl52) . This mechanism requires a non-trivial scale-dependent non-Gaussianity. 

Apart from the hypothesis of new physics, foreground contamination and instrumental systematics can break the 
statistical isotropy of the microwave background temperature map. However, these possibilities can be tested with 
the available data. Our estimates of the Cartesian components of the dipole vector, as functions of angular scale, are 
consistent for different foreground-cleaned temperature maps. The masking adopted in this analysis removes most 
contaminations from diffuse galactic emission and point sources, and our analysis procedure controls possible biases 
introduced by this procedure. In addition, realistic instrument simulations provided by the Planck collaboration 
exclude instrumental effects as the source of the observed isotropy breaking at the angular scales of interest. While 
this work was in preparation, the Planck team has made available the results of a similar analysis using the 2015 
temperature maps |58j (PLK15). Our estimates of the amplitude and direction of the dipole modulation vector on 
large-scales (/max = 60) are consistent with PLK15 analysis based on Bipolar Spherical Harmonics. The PLK15 
analysis does not provide a constraint on the scale-dependency, although shows (as for the PLK13 analysis) that the 
amplitude must decrease at higher multipoles. PLK15 shows that the amplitude of the dipole vector differs from the 
isotropic case at a level of 2cr — 3a when calculated in cumulative multipole bins [2, /max] for /max up to 320. This 
result can be compared with our geometrical test in Fig. for which similar results are found. 

Additional tests of the dipole modulation will be possible with high-sensitivity polarization maps covering significant 
portions of the sky (see, e.g., |59L 160) 1. In the standard inflationary cosmology, microwave polarization and temperature 
are expected to be only partially correlated, giving an additional independent probe of a dipole modulation; a cosmic- 
variance limited polarization map will likely double the statistical significance of the signal studied here. Gravitational 
leasing of the microwave background over large sky regions provides another nearly independent probe which will be 
realized in the near future. We will consider these possibilities elsewhere. If these probes substantially increase the 
statistical significance of the dipolar modulation signal, we will be forced into some significant modification to the 
inferred physics of the early Universe. 
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